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We use the C/N ratio as a monitor of the delivery of key ingredi¬ 
ents of life to nascent terrestrial worlds. Total elemental C and N 
contents, and their ratio, are examined for the interstellar medium, 
comets, chondritic meteorites and terrestrial planets; we include 
an updated estimate for the Bulk Silicate Earth (C/N = 49.0 ± 
9.3). Using a kinetic model of disk chemistry, and the sublima¬ 
tion/condensation temperatures of primitive molecules, we suggest 
that organic ices and macro-molecular (refractory or carbonaceous 
dust) organic material are the likely initial C and N carriers. Chemi¬ 
cal reactions in the disk can produce nebular C/N ratios of ~1—12, 
comparable to those of comets and the low end estimated for plan- 
etesimals. An increase of the C/N ratio is traced between volatile- 
rich pristine bodies and larger volatile-depleted objects subjected to 
thermal/accretional metamorphism. The C/N ratios of the dominant 
materials accreted to terrestrial planets should therefore be higher 
than those seen in carbonaceous chondrites or comets. During plan¬ 
etary formation, we explore scenarios leading to further volatile loss 
and associated C/N variations owing to core formation and atmo¬ 
spheric escape. Key processes include relative enrichment of nitro¬ 
gen in the atmosphere and preferential sequestration of carbon by the 
core. The high C/N BSE ratio therefore is best satisfied by accretion 
of thermally processed objects followed by large-scale atmospheric 
loss. These two effects must be more profound if volatile sequestra¬ 
tion in the core is effective. The stochastic nature of these processes 
hints that the surface/atmospheric abundances of biosphere-essential 
materials will likely be variable. 


terrestrial worlds | meteorites | comets | interstellar medium 

T he development of a habitable world and a stable bio¬ 
sphere requires the delivery of biogenic elements of which 
carbon and nitrogen are crucial. Carbon is the backbone for 
the chemistry of life and, in the form of CO 2 , combines with 
water to provide the greenhouse needed for a habitable Earth. 
Nitrogen is a key component of DNA, RNA, and proteins, 
while also present as the dominant constituent of our atmo¬ 
sphere. The processes that supply these crucial ingredients 
remain poorly understood. In interstellar space, C and N are 
abundant, but inherently volatile and so chiefly remain in the 
gas. Thus, the terrestrial planets, which accrete primarily 
from rocks and ices, are fed from C- and N-depleted materials 
and are carbon and nitrogen poor compared to the nebular 
disk from which they descend 1,2. The carbon and nitrogen 
depletion of rocky bodies is a general phenomenon, observable 
not just in our solar system, but in the polluted atmospheres 
of white dwarf stars, which trace the composition of disrupted 
planetesimals [3|4]. This volatile poor state of terrestrial plan¬ 
ets is partially imparted from the starting materials. However, 
further differential loss of C and N can occur due to parent 
body processes such as thermal metamorphism, core segrega¬ 
tion, planetary outgassing, and atmospheric loss [5}fTT|. 

In this work we document the evolution in the relative con¬ 
centrations of C and N from the medium (ISM) through plan¬ 
etary assembly. We show that the C/N ratio contains clues 
regarding the formation of terrestrial planets and the deliv¬ 
ery/fate of crucial volatile compounds. We first discuss the 
relative carbon and nitrogen inventories beginning with the 
ISM as it evolves to the protoplanetary disk, employing the 
perspectives of kinetic chemistry, our growing understanding 


regarding the composition of planet-forming disks, and volatile 
loss within thermally processed planetesimals. We then ex¬ 
plore the evolution of C/N during geochemical differentiation 
of a young Earth-like planet, built up by accretion of planetes¬ 
imals of a range of sizes. Key processes during this last stage 
of C/N evolution include sequestration into the metallic core 
and outgassing to the nascent atmosphere. It is of course, the 
latter that can provide ingredients for early environments and 
life, but early atmospheres are also prone to impact-driven es¬ 
cape, providing additional mechanisms for C and N processing 
and alteration. 


Carbon and Nitrogen Inventories 

Fig. [T] displays the carbon-to-nitrogen ratio observed in the 
Sun, Bulk Silicate Earth (BSE; defined in Supplementary In¬ 
formation), meteoritic samples, comets, and that inferred in 
the ISM. To provide an absolute reference frame, Fig. [^pro¬ 
vides bulk C and N abundances normalized to SiQ Key con¬ 
siderations include: 


1. In the ISM, significantly more carbon than nitrogen is in¬ 
corporated into carbonaceous dust (e.g., Fig. [21, with a 
correspondingly lower C/N ratio (<solar) in the gas. The 
primary form of condensed carbon is uncertain, but mod¬ 
els assume it is an aliphatic/aromatic hydrocarbon mix¬ 
ture fl2 13 . Interstellar ices have elevated C/N ratios, 
with a median value of C/N ~ 12 [m]. However, this does 
not account for N or N 2 , which are infrared inactive. Cor¬ 
recting for the presence of N 2 in ices (see Supplementary 
Information) leads to a C/N ratio for ices of ~1.8. 


Significance 

With the rapid pace at which exoplanets are being discovered, 
many efforts have now been dedicated to identifying which plan¬ 
ets are expected to have the ingredients necessary for the devel¬ 
opment of life. In this work we explore the relative disposition of 
the essential elements carbon and nitrogen in each stage of star 
and planet formation, using the Earth and our solar system as 
grounding data. Our results suggest that planets like the Earth 
are readily supplied with these key elements, but their relative 
amounts on the surface and in the atmosphere will be highly 
variable. 


Reserved for Publication Footnotes 


1 Further details regarding this compilation, including technique where appropriate, are found in 
Supporting Information and Table SI. 
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Fig. 1: Atomic ratios of bulk carbon relative to nitrogen in various 
solar system bodies and the ISM, based on the compilation in Ta¬ 
ble SI and references in Supplementary Information. The shaded 
ovals represent the range/errors in the measured data for each class. 
Ovals that are not closed indicate that the listed value is an up¬ 
per/lower limit. For chondrites the errors represent the range of 
determined values, while for other bodies the errors are lcr mea¬ 
surement uncertainties. For additional details an expanded version 
of this figure with labels for each point is provided as Fig. SI. 


2. Owing to low levels of observed nitrogen, particularly in 
ices (e.g., N 2 ), co mets have elevated C/N ratios when com¬ 
pared to the Sun 15 . Based on the relative disposition of 
C and N in the ISMthis can be partially achieved as comets 
accrete ISM carbonaceous dust with low nitrogen content. 
In this light, Halley is i nfer red to have more N in solid dust 
than its gaseous coma [l5j, which is also higher than esti¬ 
mates for interstellar dust (Fig. [2J) . There is also a sharp 
difference in the absolute C and N content of Halley when 
compared to the Sun-grazing comets, which have chondritic 
abundances. However, their C/N ratios are comparable. 
Anhydrous Interplanetary Dust Particles (IDP), which are 
argued to be of comet ary origin |[16 , are also characterized 
by a high carbon content, (C/SiJiDP ~ 2 17 . 

3. Carbonaceous, enstatite, and ordinary chondrites have el¬ 
evated C/N ratios compared to both the Sun and comets. 
All meteorites are carbon-poor when compared to the Sun, 
ISM dust, and comet Halley, but the key distinction is 
owing to comparative nitrogen depletion in chondrites. 
Among meteorites, carbonaceous, and enstatite chondrites 
have similar ratios, but ordinary chondrites have higher and 
more variable C/N, with mean values closer to the Bulk Sil¬ 
icate Earth. 

4. For the Earth, both the surface reservoir and Bulk Silicate 
Earth have elevated C/N ratios when compared to the Sun. 
The surface reservoir is comparable to chondritic, owing to 
preferential volcanic degassing of N 2 , leading to the N-rich 
atmosphere. There are two other estimates of the C/N ra¬ 
tio for the Earth, C/N ~ 53=L20 [l and 425=L329 [7] (see 
Supplementary Information); in each case the estimated 
C/N ratio of the BSE is quite high[^] 

These points tell a tale of an increasing C/N ratio from 
the cold ISM, to cold cometary formation zones, to the inner 
nebula to planets. While there is a broad overall trend, it is 
also clear that local processes can lead to variations within 
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Fig. 2: (a) Bulk carbon-to-silicon and (b) nitrogen-to-silicon ra¬ 
tios. Symbols are identical as in Fig. 1. References and an expanded 
figure with labels for each point in Supplementary Information. 


and between certain classes. Whilst much of this secular vari¬ 
ation must be due to high N volatility, additional processes 
including preferential dissolution in molten Fe,Ni-metal may 
contribute during assembly of planetesimals and planets. No¬ 
tably the, BSE has the highest C/N ratios when compared 
to most solar system material, except for ordinary chondrites. 
This high C/N ratio, coupled to overall depletion of the BSE 
in absolute C and N concentrations, indicate that processing 
and loss of these key volatiles occurred during planetary as¬ 
sembly. Thus, there are several processes active early in the 
young solar system that control how planets are seeded with 
C and N. Below we explore the processes responsible for these 
trends, which must relate to the supply and fate of two of the 
key ingredients of our biosphere. 

C/N - ISM to Solar System Rocks and Ices 
Elemental Condensation/Sublimation. The depletion of car¬ 
bon and nitrogen within meteorites and Sun-grazing comets 
is partially related to the high volatility of primary carriers 
within ices (e.g., CO, N 2 ), as opposed to the carbonaceous 
dust component. In the case of C, it is further complicated 
by the fact that interstellar carbon grains provided to the disk 
need to be destroyed at some level. This must be the case as 
the C/Si ratio of meteorites (and Sun-grazing comets) is an 
order of magnitude below that of the ISM dust (e.g., Fig. [2|. 
Carbonaceous dust destruction mechanisms likely include re¬ 
actions with atomic O and OH (oxidation), photodissociation 


2 Limited estimates for Mars and Venus are supplied in the Supplementary Information. 

3 In the supplemental material we discuss the available information regarding the formation loca¬ 

tion^) of Halley and the Sun-grazing comets. 


2 | www.pnas.org/cgi/doi/10.1073/pnas.0709640104 
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Fig. 3: Fraction of elemental C (left panel) and N (middle panel) in molecular gases or ices from a kinetic chemical mode (age=l Million 
years), (right panel) C/N ratio predicted for organic ices, the total ice mantle (organics and other volatiles shown in the other panels), and 
for refractory carbon (carbonaceous dust). For C the evolution is complicated by the uncertain fate of carbonaceous dust which carries 
~50% of elemental inventory. Based on the composition of meteorites, some dest ruction of carbon-rich dust is required. The red shaded 
area represents the (uncertain) radii where carbon grains might be destroyed [19 . Within this zone we artificially reduce C incorporated 
in carbonaceous dust to zero. The principal product of the destruction is CO, "which is artificially shown in the figure as CO 7 . 


of small (< 50 C,H atom) clusters [l8]— 20 , and parent body 
processing (see following section). Oxidation can destroy car¬ 
bon grains, but not silicates, but requires high gas temper¬ 
atures (>few hundred K). This process can be active on the 
surface of the inner disk and/or during the hot, early evolution 
stages leading to C depletion in solar system rocks P1P1- At 
greater distances from the Sun, or at later times, interstellar 
carbon grains need to be significantly preserved, or perhaps 
gain additional carbon by extraction from CO [21], to match 
the carbon content of comet Halley and IDPs. The difference 
between comet Halley, IDPs, and the Sun-grazing comets may 
relate to the processes discussed abovej^j 

To explore this chemical evolution, we present the main car¬ 
riers of elemental C and N within a generic protoplanetary disk 
kinetic chemical model (see Supplementary Information) in 
Fig.0 The temperature of microscopic solids from 0.5-3 AU 
will be uniformly above the CO condensation temperature of 
~20 K. Nitrogen, primarily N 2 , is also gaseous here, though 
NH 3 (condensation temperature of ^80 K) has lower volatility 
than the main carbon ice carriers (CO, CO 2 ) and could pro¬ 
vide nitrogen to some solids. In this example, it is only inside 
2 AU that NH 3 ice begins to sublimate. Organics are present 
at these radii and, due to gas-phase chemical processing some 
of the C in CO and N in N 2 can be reprocessed in the gas 
leading to organic ices [2T . 

Thus, as seen in Fig. [3] the carbon carriers with the lowest 
volatility and the primary carriers of C in inner solar system 
solids will be organic ices and/or carbonaceous dust grains. 
These are the carriers of elevated C/N ratios into the warmer 
inner disk. We provide a brief discussion of the relative sub¬ 
limation temperature of organic materials in the Supplemen¬ 
tary Information. In the comet-forming zones (R > 5 AU), 
carbon is primarily found as carbonaceous dust and as CO ice 
(>20 AU) or CO vapor (<20 AU). At least 50% of the carbon 
resides in dust with the condensation/sublimation tempera¬ 
ture > 400 K. For the case of nitrogen, in this model the inner 
10 AU are dominated by NH 3 ice and N 2 gas, while beyond 
10 AU there is even partitioning between N 2 ice and NH 3 ice. 

Based on this analysis, the C/N ratio in the starting mate¬ 
rials for planets depends on the mixture of carbon in carbona¬ 
ceous dust and organic ices coating dust grains. Given the 
low carbon content of chondritic meteorites, the presence of 
interstellar carbonaceous dust grains might be time/distance 


variable in the disk. In contrast, organic ices might be created 
in situ at all radii with C/N of 1-6 (Fig. [31, while incorporat¬ 
ing some carbonaceous dust would give higher values. Comet 
Halley provides a limit for the carbonaceous dust and ice con¬ 
tent in the ice giant planet-forming zone, with C/N = 3.7-11.5. 
Disk chemistry can thus implant C/N ratios ranging from 1 
to ~12, values consistent with comets and the low end of the 
primitive chondrites. 

Loss During Parent Body Processing in the inner solar sys¬ 
tem. A key issue in interpreting the C/N ratios of meteorites 
is that they are not pristine records of the dust and ice that 
accreted to form planetesimals in the solar nebula. Meteorites 
may provide a good analog to such processes, as their parent 
bodies have been heated and differentiated, modifying both 
their C and N concentrations. |5|6| . Chondritic meteorites, are 
believed to be direct aggregates of dust grains from the solar 
nebula, but were exposed to high temperatures (up to 1300 K) 
and experienced different degrees of thermal metamorphism, 
in some cases resulting in elevated C/N. This does not mean 
that high C/N meteorites, such as ordinary chondrites (Fig. 1) 
were direct sources of high C/N materials delivered to Earth, 
as chemical and isotopic argu men ts rule out these samples as 
primary terrestrial precursors [22]. Rather, the processes that 
they experienced may have been analogous to processes that 
occurred on parent bodies of objects from which the Earth 
accreted. 

Volatile species, including some organics, would decompose 
in meteorite parent bodies at high temperatures and the lib¬ 
erated C and N could migrate to the surface as a gas and be 
lost to space, or could be incorporated into metallic Fe,Ni (ka- 
macite and taenite) owing to their siderophilic behavior [23] . 
Thus the abundances of both C and N are expected to diminish 
within bodies that were exposed to higher temperatures and 
more prolonged heating. Thi s is seen in both ordinary 24-26 
and carbonaceous chondrites [27]; additional evidence isTound 
in the volatility of organic components and their carbon and 
nitrogen isotopic ratios [28] . 

The extent of C and NToss during parent body processing 
is determined by the initial compounds accreted, the chem¬ 
istry that occurs during heating, and the intensity of heating. 
For example, the vigor of C and N release from organics dur¬ 
ing thermal metamorphism depends on the redox conditions, 
particularly for C, which is incorporated into a gas under oxi- 
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dizing conditions, but is s tabi lized in refractory compounds at 
more reducing conditions [23]. This could account for much of 
the variation in C/N ratio seen in ordinary chondrites. How¬ 
ever, enstatite chondrites have most of their N in the form 
of nitrides while carbon is found chiefly as graphite and car¬ 
bides [29]. The reduced character of enstatite chondrites likely 
stabilizes solid C and N compounds to higher temperatures, 
inhibiting large-scale modification of the original C/N value. 

Continued heating on parent bodies would lead to melting, 
and the formation of iron cores surrounded by silicate mantles 
and any C and N not lost to space would be distributed accord¬ 
ing to their low-pressure metal/silicate partition coefficients. 
Growing planets accrete both silicate and metal portions of 
such a body, either in individual events or, if fragmentation 
occurs 30 , in aggregate of multiple accretion events. How¬ 
ever, it is chiefly the volatile loss that occurs during heating 
that would determine the final inventory of materials delivered 
to a growing planet. Thermally processed ordinary chondrites 
may thus show us the type of processing that occurred on 
the planetesimals and embryos accreted to the Earth: ther¬ 
mal evolution of small bodies derived from nebular dust and 
ice produced differentiated bodies with enhanced C/N ratios. 


C/N - Forming Terrestrial Worlds 

Volatiles accreted to the Earth influence climate and biosphere 
development if they are available to the near-surface environ¬ 
ment, either as remnants of the early primordial atmosphere 
or from material incorporated into the silicate Earth and later 
out gassed via volcanism. Because the early atmosphere was 
susceptible to loss via impacts Em the latter may be most 
important. But, volatiles containeam silicate could be drawn 
down into the core during the era of silicate-metal segregation. 
Consequently, the volatiles presently in the BSE have escaped 
loss to space or sequestration by the core, leaving them poten¬ 
tially available to the modern surface environment. 

Current models of early planetary differentiation postulate 
that formation of the core and of an early massive atmosphere 
are associated with a deep magma ocean, with metal-silicate 
partit ioning established at high mean pressures and tempera¬ 
tures [33||34 - possibly at conditions that became more oxi¬ 
dizing with time 35||36 , and with a significant fraction of the 
planetary volatile budget released from the molten silicates 
into the atmosphere [37[[38] . Both C and N are moderately to 
highly soluble in molten Fe-ri ch alloy (819,39 and compara¬ 
tively insoluble in silicate [40||42]. Thus, much of the Earth’s 
initial C and N inventories likely segregated to its metallic core 
or were released to the atmosphere. The quantitative effects 
of these processes on BSE C and N content, and on the C/N 
ratio, depends on partitioning between molten Fe-rich alloy 
and molten silicate, which in turn depends on the conditions 
of metal-silicate exchange and on the solubilities of C and N in 
silicate at the magma ocean-protoatmosphere interface. Also, 
the effects of core segregation depend on the fraction of metal 
equilibrated with volatile-bearing silicate [43], as much of the 
core may have segregated prior to bulk volatile delivery |44| . 

Experimental studies establish that C and N partition 
strongly into molten Fe-rich alloy coexisting with silicate melt. 
Though values for both elements vary with experimental con¬ 
ditions and compositions, C partitions more strongly into 
metal than does N. Characteristic mass partition coefficients 
for C (Dj? et / S * z ) are 10 2 -10 4 [8 10 45 47 , compared to those 
for N (D^ et ^ sd ) of ~20 [ 9 ]. The solubilities of C and N de¬ 
pend strongly on oxygen tugacity, with C more soluble than 
N in silicate melt under oxidizing cond itions, and the opposite 
under highly reducing conditions . . _ 

Although terrestrial volatile accretion and differentiation 
were dynamic processes evolving over a range of conditions, 



Fig. 4: Calculated C/N for the BSE resulting from a whole-mantle 
magma ocean equilibrated with an overlying atmosphere and vari¬ 
able amounts of core-forming metal |43| . All calculations assume 
the proto-Earth begins with C/N = 25 and these elements are ini¬ 
tially present in the silicate. (A) shows C/N for a BSE deriving 
from the sum of the magma ocean and the atmosphere, whilst (B) 
shows the BSE C/N derived only from the magma ocean, assuming 
that the primordial atmosphere is lost to space. (Oxidized): C 
in vapor is assumed to be chiefly CO 2 and N is dissol ved as N2. 
(Reduced): C in vapor is assumed to be chiefly CO 47 . N is 
dissolved as N 2 . (Very Reduced): C is highly insoluble and N 
highly soluble. Quantitative solubilities for case 3 are poorly con¬ 
strained, but for C and N must be very low and high, respectively. 
See Supplementary Information for calculation details. 


insight can be gained from simple models that begin with 
a molten mantle in equilibrium with some fraction of core¬ 
forming metal and with an overlying atmosphere |37[|43| . To 
track the evolution of the C/N ratio during planetary process¬ 
ing, in Fig. [4] we consider two end-member cases, both adopt¬ 
ing an initial stage of chemical equilibrium among a fraction of 
core-forming iron-rich alloy or metal, a molten silicate mantle, 
and an overlying atmosphere, followed by the segregation and 
isolation of the core from the mantle. In case A, the atmo¬ 
sphere returns to the silicate mantle to form the BSE, whereas 
in the other, case B, the atmosphere is blown off to space and 
lost from the Earth. Motivated by the earlier discussion, we 
assume an initial C/N ratio of 25. 

The most oxidizing condition (orange curves, Oxidized) is 
a scenario in which carbon in vapor is assumed to be CO 2 
and N in vapor as N 2 . We first discuss Case A where the 
magma-silicate Earth is enshrouded by a primordial atmo¬ 
sphere. Here, the initial equilibrium leaves a large portion of 
N in the atmosphere with most C in the condensed phases and 
predominately in the core. Thus, if the primordial atmosphere 
was retained, the BSE will end up with a C/N ratio slightly re¬ 
duced from the initial value, provided a small fraction of metal 
equilibrated with the mantle. The ratio approaches zero if the 
entire core (metal/silicate~0.5) equilibrated with the molten 
silicate and sequestered most C away. In case B (atmosphere 
blow off) the BSE C/N ratio would be elevated with respect 
to the initial value if core sequestration is limited, but reduced 
if a large fraction of the core equilibrated with the mantle and 
strongly depleted the C of the molten silicate. 

The second scenario (blue curves, Reduced) represents mod¬ 
erately reducing conditions wherein carbon in vapor is as¬ 
sumed to be CO and nitrogen as N 2 . If the atmosphere is 
retained, the results are qualitatively similar to the first case, 


4 | www.pnas.org/cgi/doi/10.1073/pnas.0709640104 
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but if the atmosphere is blown off this results in a BSE C/N 
ratio that does not exceed the initial value even if no C was 
sequestered by the core. The lower than initial C/N ratio is 
derived from the smaller solubility of C in molten silicate. 

The final scenario (black curves, Very Reduced) represents 
a highly reducing condition where C is assumed to be highly 
insoluble, with a speculated solubility of 0.05 ppm/MPa, and 
N highly soluble, with a speculated value of 50 ppm/MPa. 
Contrary to the previous two cases, initial equilibrium among 
metal, silicate, and atmosphere leaves most N in the condensed 
phase and a significant portion in the silicate, whereas C re¬ 
sides primarily in the atmosphere. For case A, the BSE C/N 
ratio is slightly elevated from the initial value because some N 
was sequestered by the core and C escaped core sequestration 
by being part of the primordial atmosphere. On the other 
hand, for atmospheric blow off (case B), most C would be lost 
and the C/N ratio of the BSE would be markedly reduced. 

All models assume that £)™ et / s ^ jj™et/sii gQ conse _ 
quences for the BSE C/N are affected by: (a) the relative sol¬ 
ubilities of C and N in silicate melt, and (b) the proportion of 
metal that equilibrates with the magma ocean. As illustrated 
by Fig. 13 if C is modestly soluble in silicate and N insoluble, as 
is true under oxidizing or modestly reducing conditions, then 
C is chiefly dissolved in condensed phases, and predominantly 
so in metal. The BSE (silicate + atmosphere) thus has greatly 
diminished C/N, even for small fractions of available molten 
alloy. If, however, the amount of metal is small, then following 
atmospheric loss, the residual silicate mantle may have high 
C/N. Alternatively, if C is highly insoluble and N very solu¬ 
ble, as may be true under very reduced conditions 41 47 , then 
the proportion of N dissolved in condensed phases, including 
core-forming metal, is enhanced. The atmosphere is therefore 
carbon-rich, and so the BSE can have a high C/N ratio. 

In sum, under most plausible scenarios, magma ocean- 
related core formation, combined with atmospheric blow off, 
produces much lower C/N ratios than observed. The excep¬ 
tions is a speculative circumstance of highly reduced condi¬ 
tions where N is more soluble than C in silicate. More likely, 
N-rich atmospheric blow-off is necessary to elevate the C/N 
ratio from carbon depleted levels (Fig. 0) to the current value. 

Discussion 

The C/N ratio provides diagnostic information regarding 
physical processes active during the earliest stage of planetary 
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Supporting Information 

Bulk C and N Abundance Compilation and Estimation 

Below we summarize the information compiled and estimated 
for Figures 1 and 2 in the main manuscript. For completen ess 
we provide expanded versions of these figures here in Fig. [Si] 
and Fig |S2| Table SI provides our compilation in tabular form. 


Terrestrial Worlds. The Bulk Silicate Earth (BSE) assumes, 
for simplicity, that the Earth is comprised of two components: 
the core and the silicate mantle. In this form the BSE is the 
entire Earth, including its fluid envelopes (oceans and atmo¬ 
sphere) minus the core. Previous estimates of the C/N of the 
BSE differ markedly. Marty [7] estimated a C/N molar ratio 
of 365±282, whilst in contrast Halliday’s preferred model 1 
amounts to 46±17. This factor of 8 difference led us to re¬ 
visit constraints on BSE C/N, using an approach somewhat 
different from previous work. 

The BSE C and N inventories derive from a combination of 
the surface and mantle reservoirs (Table SI). For the surface 
reservoir, we take total carbon to be 9.3 ±0.9 xlO 22 gm 


and total nitrogen to be 6.4±1.1 x 10 21 gm 50 , yielding a C/N 


49 


molar ratio of 16.9±0.3. For the interior reservoir we consider 
two separate sources, based on the well-established assumption 
that the mantle consists of a gas-depleted source (DM) that 
produces mid-ocean ridge basalts and a gas-enriched source 
(EM) related to oceanic island basalts [51 . 

We assume that the deple ted reservoir constitutes between 
30-80% of the mantle mass [52], and estimate the C concen¬ 
trations of each reservoir based on C/Nb and C/Ba ratios of 
oceanic basalts compiled by Rosenthal et al. (in press), yield¬ 
ing 20±7 and 165±55 ppm for the DM and EM, respectively. 

Marty \7 and Halliday [T| estimated mantle N concentra¬ 
tions in part from the esti mate of mantle C/N of 535±224 
Marty and Zimmerman [53]. However, this ratio derives from 
gas liberated from basalt crushing experiments. Because N 2 
fractionates from CO 2 substantially during vesiculation, C/N 
ratios of basalts and their vapor bubbles vary by more than 
two orders of magnitude [53,54 , complicating reconstruction 
of source C/N. Therefore, we estimate the N concentrations 
of each reservoir based on N/Ar ratios. Unlike C/N, N/Ar 
ratios of vesicles in basalt do not change significantly owing 
The mass of Ar in the mantle is con- 
inventory of K, 280±60 ppm, which over 

. Sub- 
of 6.6 


to degassing 153,54] 
strained by the BSE 

4.55 Ga has produced 1.55±0.33 x 10 2U gm of 4U Ar 
tracting atmospheric |56 and crustal 57 Ar masses 


55 


and 0.9 xlO 19 gm, respectively, leaves 8.00 ±3.3 X 10 iy gm of 
Ar in the present-day mantle 55 . Because the N/Ar ratios 
of the depleted (DM) and enriched mantle (EM) reservoirs 
are see ming ly distinct (N/Ar = 74 ± 56 in DM and 105±35 
in EM; [58]), conversion of the bulk mantle Ar mass to the 
bulk mantle N requires estimates Ar concentrations in each 
reservoir. Taking the DM to have 3.3±1.6 ppb Ar (Marty and 
Dauphas [58], but also consistent with Allegre et al. [5H and 
Ballentine [59] ) allows calculation of the concentration of Ar 
in EM by subtraction from the bulk mantle Ar reservoir. The 
above information is sufficient to calculate the C/N ratios of 
the principle terrestrial reservoirs (surface reservoir, bulk man¬ 
tle, BSE) (Table SI), taking in to account the uncertainties 
in individual concentrations and ratios must be that one ac¬ 
counts for the uncertainties from each constraint. T his i s done 
with a Monte Carlo simulation, as illustrated in Fig. |S3| which 
yields C/N ratios for the bulk mantle and BSE of 80 and 46±9, 
respectively. The mantle ratio is reported as a lognormal dis¬ 
tribution because the arithmetic mean and standard deviation 
are strongly affected by simulations that produce very small 


mantle N (and thus extreme C/N ratios). This BSE C/N 
estimate is in nearly exact agreement with the previous esti¬ 
mate by Halliday H, 46±17, and markedly lower than that by 
Marty [ 7 ]. 

Elemental abundances for the atmospheres of Venus and 
Mars are taken from the compilation of Hailliday et al. 1 . 
The atmospheres of Mars and Venus have C/N comparable 
to the Earth’s surface (Supplementary Fig. SI), though in¬ 
ferences for these planets is limited by lack of constraints on 
interior reservoirs. In the small sample of Martian meteorites 
the C/N ratio appears to vary, but are elevated above the solar 
value [29] . 


Meteorites. Chondritic meteorites are a class of meteorites 
that escaped differentiation. Their bulk compositions gener¬ 
ally match that of the solar photosphere, except for the most 
volatile elements. These meteorites are broken into groups 
based on their bulk chemical, mineralogical, and oxygen iso¬ 
topic compositions. In addition, the proportion and sizes of 
the different nebular solids they contain can be used to dis¬ 
tinguish among the groups. These groups are also categorized 
into three distinct classes of chondrites: the Enstatite, Or¬ 
dinary, and Carbonaceous chondrites (EC, OC, and CC, re¬ 
spectively). While chondrites are thought to be aggregates 
of solar nebula dust, these meteorites have experienced dif¬ 
ferent degrees of planetary processing in the form of thermal 
metamorphism and aqueous alteration. Aqueous alteration 
was most pervasive in the Cl and CM chondrites, resulting 
from the melting of ice and reaction of rock with liquid wa¬ 
ter at temperatures <400-500 K. While evidence for the ac¬ 
tion of water in other chondrites has been reported, thermal 
metamorphism, with meteorites reaching temperatures of up 
to 1300 K, is much more common in the EC, OC, and some 
of the other CC (notably CV and CO) meteorites. 

For meteoritic C/N ratios we have used the values deter¬ 
mined by Alexander et al. J60l[6l] for Cl, CM, CR, and CO 
chondrites, while the data for CV chondrites are taken from 
the work of Pearson et al. [27]. For these data points we have 
determined the weighted mean and errors based upon the in¬ 
dividual meteoritic measurements. Thus the error represents 
the range of the data relative to the errors of the individual 
measurements. For carbon and nitrogen elemental ratios rela¬ 
tive to silicon we have used the bulk abundance measurements 
of Wasson & Kallemeyn 162 , with an assumed uncertainty of 
20% which is based upon their discussion of the relative un¬ 
certainties within carbonaceous chondrite groups. We also ex¬ 
tended the error bar range based on the variation seen in the C 
and N content within a given class from the references above. 
For ordinary chondrites, information was gathered from differ¬ 
ent sourc es often studying the same meteorite but in different 
elements 24] 25 63 65 . The error reflects the range of es¬ 
timated C/N ratios for each class. Abundances relative to 
silicon were estimated using the range of C and N from the 
previous sources for a given class and the average Si abun¬ 
dance for that class taken from Fodders V Fegley [66]. For 
both ordinary and carbonaceous chondrites the error bars thus 
represent the variation of C and N rather than that of Si. 


Comets. Abundances for comet Halley rel y o n compila¬ 
tions/analysis of Fomenkova 67 , Jessberger 68], Delsemme 
69 , and Wyckoff et al. [l5]. In situ analysis ofthe refractory 
carbonaceous component was performed by instruments on 
board Vega-1, Vega-2, and Giotto spacecraft that performed 
fly-by’s within the coma and that, with further analysis, pro¬ 
vided elemental abundances. For carbon, nitrogen, and sili¬ 
con, Fomenkova [67] state that gas and dust abundances of 
Halley for carbon are solar, nitrogen is a factor of two be- 
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low solar, while silicon is within a factor of 2 of solar. Jess- 
berger [68] states that the rock forming elements are within a 
factor oT two of solar and chondritic; thus we adopt a factor 
of 2 error on the Si abundance which then dominates that cal¬ 
culation. For carb on we assume that 36% is bound in molec¬ 
ular volatiles 69], while only 10% of the nitrogen inventory 


is tied up in such species 


niy 11 

i, 


Assuming solar abundances 


the C/N (gas and dust) ratio is (C/N) 0 /2 and we adopt a 
50% error. The evaporation of Sun-grazing comets provide 
another unique estimate of bulk cometary abundances. Two 
comets have measurements of C and Si (C-2003K7 and C/2011 
W3:Lovejoy) [70 71 , while Comet Lovejoy has an additional 


measurement ornitrogen (provided by J. Raymond and given 
here). Errors are assigned a value of 50% which is based upon 
the estimated range in the C/Si ratio for C-2003K7 [70] and 
we extend this to the other measurements. 

Regarding the formation zone of Halley and the Sun-grazing 
comets, we note that the two Sun -gra zing comets could be 
fragments of the same progenitor [72]. However it is likely 
that the larger parent (s) originated within the Oort cloud 73] , 
which itself is populated from comets formed near the gi¬ 
ant planets and perturbed to larger orbits 74 . The origin 
of Halley-type comets has been a matter of debate; however 
the recent detection of a trans-Neptunian object in a highly 
inclined retrograde orbit suggests a source region for Halley- 
type com ets with a primordial origin at large distances (> 50 
AU) [75. 


Interstellar Medium. In the ISM, elemental abundance mea¬ 
surements can be obtained via spectroscopic techniques to¬ 
wards low density (<1000 cm -3 ) atomic gas in front of bright 
stars. These gas-phase data are compared to some refer¬ 
ence abundance measurements, either the Sun or young stars 
(which likely represent the original composition of gas in their 
galactic location). For a given element the deficit in abundance 
from the gas when compared to the solar/stellar standard is 
assumed to be sequestered into dust grains. The presence of 
dust along these lines of sight is confirmed by its differential 
absorption effects on the broad band stellar spectrum. 

Thus, Jensen et al. [76] surveyed 30 sight lines in transitions 
of atomic nitrogen and compiled abundances from previous 
work. They find a gas phase nitrogen abundance at high ex¬ 
tinction (Ay > l m ) of 49 ± 4 ppm (relative to H) and, using 
adopting B stars as the standard, find that 14 ppm is missing 
from the gas phase. This missing nitrogen could potentially 
be incorporated into the carbonaceous dust component and 
represents our limit for nitrogen in the solid state. 

Based on elemental abundances measured toward sight 
lines that intersect diffuse interstellar gas/dust, the wave¬ 
length dependent extinction of starlight, and spectral fea¬ 
tures with unidentified carriers (e.g. the 2200 A extinction 
bump and mid-IR emission bands), it is a thought that sig¬ 
nificant percentage of interstellar carbon resides in some re¬ 
fractory component. In 5 sight lines with hi gh a verage density 
(logio (nH/cm -3 ) > 0.5), Parvathi et al. 
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estimate that 

371 ± 37 ppm (relative to H) of carbon resides in interstellar 
dust grains, assuming B star standard abundances. Along the 
same, higher density sight lines the abundance of gas-phase 
carbon is 94 ± 13 ppm. The recent dust models of Jones et 
al. [l2] and Chiar et al. 13 adopt 233 ppm and 150 ppm, 
respectively. For our purposes it is clear that there is a large 
amount of carbon in interstellar dust grains and we adopt 250 
± 50 ppm to account for the disparate abundance estimates. 

As noted in the main te xt t he median value of C/N ratio 
of interstellar ices is ~ 12 [ 14]. These measurements do not 
account for N or N 2 , which are infrared inactive. Based on the 


similar sublimation temperature between CO and N 2 78 , we 
assume (CO/N 2 )i ce ~ 1. Further, assuming all N is inNF and 
50% of C in CO, gives a C/N ratio for ices of ~1.8. 


Assumptions for Equilibrium Calculations 

In our calculations in Fig. [4] we made the following assump¬ 
tions. Oxidized: — 500 , jj™ et / sxl — 20, C in vapor is 

assumed to be chiefly CO 2 , with d epth -aver aged C solubility 
in silicate liquid of 1.6 ppm/MPa 40, N is dissolved as N2, 
with depth-averaged solubility=l ppm/MPa |41|. Reduced: 
jjrnet/sii _ ^qqq jjfriet/sii _ q j n va p 0r [ s assumed_to be 

chiefly CO, with solubility equal to 1 ppm CO/MPa [47]. N 
is dissolved as N 2 , with solubility=l ppm/MPa [4l]. Very 
Reduced jj™ et / sd — 1000, jj 1 ^ et / s ' 11 — 20, C is highly insol¬ 
uble (0.05 ppm/MPa) and N highly soluble (50 ppm/MPa). 
We note that quantitative solubilities for case 3 in the figure 
lack constraints; for C and N they must be very low and high, 
respectively. 


Kinetic Chemical Model 

Fodders [79] explored the chemistry in the solar system within 
the framework of the 50% condensation temperature for ther¬ 
modynamic equilibrium. This concept implies that the chemi¬ 
cal species produced in the disk are independent of the starting 
materials. However,as discussed by Lodders, the history and 
carriers matter for reactive volatile elements. Thus, the 50% 
condensation temperature is not a useful guide in this case, as 
thermodynamic equilibrium will not be achieved and we adopt 
a kinetic model to explore the evolution of these elements. 

The kinetic chemical model presented in Fig. 3 of the main 
text is a compilation of ch emi cal abundance predictions from 
the model of Cleeves et al. [80 . The physical structure is a 2D 
azimuthally symmetric model that is motivated by the obser¬ 
vations of solar analog protoplanet ary disks. Greater details 
regarding the exact physical prescription can be found in the 
references listed above. For completeness the basic physical 
structure is provided below in Fig. |S4| Our methodology is to 
take an adopted dust density structure constrained by obser¬ 
vations. The gas density in this model, when compared to the 
dust, has a vertically integrated gas to dust mass ratio of 100. 
However, in practice the gas to dust mass ratio is higher in the 
upper layers and reduced in the midplane to mimic the effects 
of dust settling. With the gas and dust density prescribed we 
determine the dust temperature structure using the radiative 
equilibrium code TORUS [81]. For this calculation the central 
star is assumed to have mass of M = 1.06 M 0 , R = 1.83 R 0 , 
and an effective temperature of 4300 K. Due to high densities 
and frequent dust-gas collisions, the gas temperature is cou¬ 
pled to the dust for much of the disk mass, but is elevated at 
the surface due to stellar X-ray heating. Within the frame¬ 
work provided by the dust and gas density we determine the 
UV and X-ray photon fields. For the input spectra in the UV 
we use the spectrum measured towards TW Hya, the closest 
star/disk system. The X-ray spectrum is estimated assuming 
a two-temperature optically thin plasma with one component 
at T = 9 MK and another at 30 MK, and propagated as¬ 
suming X-ray photoabsorption cross-sections. The resulting 
distribution of UV and X-ray flux is then calculated using the 
Monte-Carlo model developed by Bethell & Bergin [82]. For 
more explicit details and complete references the reader is re¬ 
ferred to Cleeves et al. [80] , 

The chemical model includes gas-phase processes, gas-grain 
interactions, and a limited set of grain surface reactions for 
the main carriers of C, O, and N linked through a network 
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of ~5900 reactions assuming initial chemical abundances 80 . 
To create Fig. 3 we ran the model for a million years and com¬ 
piled the predicted abundances. Based on these predictions we 
summed up the abundances of all carriers of carbon and nitro¬ 
gen to isolate the main carriers. These are (both gas and ice 
phase): CO, CO 2 , N 2 , NH 3 , and organics. Refractory carbon 
contains ~50% of elemental carbon in the interstellar medium 
and we assume it is oxidized [l 9 interior to ~2 AU where it is 
released as CO (listed as CO' in Fig. 3). The location in the 
disk where this process is active is highly uncertain. The main 
organic carriers of C, where the abundances were summed for 
the figure (separately as gas and ice) are: HC 3 N, C 3 , HCN, 
CH 3 OH, C 3 H, HNC, C 4 H 2 , C 2 H 4 , OCN, C 3 0, H 2 CO, and 
CH 4 ; for nitrogen the organic carriers are: HC 3 N, HCN, HNC, 
and H 2 CN. An important issue is the presence of NH 3 in the 
ices even at 2 AU. This is due to the relatively high sublima¬ 
tion temperature of NH 3 compared to N or N 2 (at relevant 
pressures, NH 3 : ~ 80 K for pure ice or higher if bound with 


water; N 2 : ~ 20 K) [78|8 3 . In this light the depletion of nitro¬ 
gen in comet ary ices is difficult to interpret. However, a recent 
exploration of nitrogen in protoplanetary disks can match the 
nitrogen content in cometary ices provided nitrogen is initially 
supplied to the disk as atomic N, which is then converted to 
gaseous N 2 [84] . 

The estimated sublimation temperature of organic ices de¬ 
pends on the binding stre ngth to the surface; based on lab¬ 
oratory measurements [85], this is on the order of 150-250 K 
for low molecular weight species (the latter corresponds to a 
binding energy of 7500 K, for a molecular mass of 60 amu, and 
a pressure of ~10 -4 bar). For carbon in organic “refractory” 
form we can use the literature on sublimation enthalpies for 
hydrocarbons which are ~ 100 kJ/mol [ 86 ], with a correspond¬ 
ing sublimation temperature of ~ 400 K. Some fraction of the 
C may be present as macromolecular, cross linked species, and 
would survive to even higher temperatures. 
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Fig. SI: Expanded version of Fig. 1 in the main manuscript. Ratio of bulk (or as noted) carbon relative to nitrogen in various solar 
system bodies and ISM. Errors denoted with a diamond represent ranges within the samples and not the measurement error. Normal 
errorbars represent the lcr uncertainty in the measurement. References for this compilation are given here. 
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Fig. S2: Expanded version of Fig. 2 in the main manuscript, (a) Bulk elemental carbon to silicon atomic ratio and (b) nitrogen to silicon 
atomic ratio. Errors denoted with a diamond represent ranges within the samples and not the measurement error. Normal errorbars 
represent the lcr uncertainty in the measurement. References for this compilation are given here. 
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Fig. S3: Results of Monte Carlo simulation showing compositions of surface reservoir, bulk mantle, and bulk silicate Earth 
(BSE). Concentrations are in ppm by weight, calculated from the mass of C and N in each reservoir divided by the mass of the 
BSE (4 xlO 27 gm). The Monte Carlo simulation was performed with 1000 individual calculations. 
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Fig. S4: Figure illustrating the adopted disk physical structure for the kinetic model as a function of radius (R) and height 
(Z). This includes, the gas density (a), dust density (b), d_ust temperature (c), gas temperature (d), UV photon flux (e), and 
X-ray photon flux (f). Figure taken from Cleeves et al. 
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Table SI: Compilation of C and N Content in Solar System Bodies 
Relative Atomic Ratio* 

Source C/Si N/Si C/N Notes on CN Ratio 

Terrestrial Planets 

Earth (BSE) l.ips)!^!^ 2.5(-5)+°[l®f 49.0 ±9.3 This Text* 

Earth Surface 2 - 5 (- 4 )- 2 l(- 5 ) 14.3 ±2.9 This Text* 

Earth Mantle 8.7(-4)t^[l^ 9.2(-6)±4;®|l®j 80 ±35 This Text* 


Venus Atm. 

3.8(—4) 

2.8(—5) 

13.5 ±2 ] 

L t 

Mars Atm. 

1.6(-7) 

8.8(—9) 

18.0 ±2 ] 

L f 


Meteorites* 


Cl 

0 71 + 0 - 14 

U * 1 1 -0.50 

o 

o 

4^ 

to 

1 + 
oo 

oo 

17.0 ±3.0 

Cll 

60 


CO 

o 

o 

-'T 

1 + 
oo 
oU 
to 

O 

O 

o 

CO 

1 + 
oo 

oo 
oo 
to to 

12.2 ±3.2 

C03 61 


CM 

o.4o±° : °f 

0.0221°;°}° 

21.7 ± 1.3 

CM2 6( 

D 

cv 

O 

O 

OO 

1 + 
oo 

OH 

0.0041°;°°! 

20.8 ±3.7 

CV3 27 


L 

O 

O 

-4 

1 + 
oo 

oo 

0.00ll°;o°f 

45.2 ±25.2 

L3 

24 : 

>5 64 65 

LL 

O 

O 

-4 

1 + 
oo 

oo 

0.00ll°;°°f 

56.0 ±43.7 

LL3.5 k 


H 

oo 
do 
+ 1 

O 

O 

0.00013I'q 0 oooos 

452.0 ±298.0 

H3,4 \2 l 

ppt 

EH 

0 f)9?+ 0 - 008 

U.UZZ_Q 008 

O 

O 

o 

to 

1 + 
oo 

oo 

oo 

13.7 ± 12.1 

EH5.6 : 

29 


Comets 


Comet Halley (dust) 

5.7+0 

1 0+1.0 

4 ’ U -0.5 


15 

67 

69 


Comet Halley (dust ± gas) 

8.3±|;1 

1 1+ 2 - 1 

4 ' 4 — 0.5 

7.5±3.75 § 

67 




Comet C/2011 W3 

0.2ll°;}° 

O 

o 

CO 

l + 
oo 

oo 

6.2±0.9 

71 

If 



Comet C/2003 K7 

O 

o 
00 
l + 
oo 

oo 



71 




Interstellar Medium 

ISM (Gas) 

41.71!;° 

21-9111 

1.9±0.3 

76 

77 



ISM (Dust) 

5-851};! 

< 0.33 

> 17.57 

12 

13 

76 

77 



*See text for discussion of the upper and lower bounds, which represent ranges in measurements and/or uncertainties in a compiled 
measurement. 


11.0( — 4) corresponds to 1.0 X 10 4 

£ As noted in the discussion regarding meteorites, the C/N ratio in carbonaceous chondrites was determined from different references than 
the C/Si and N/Si. However the measurements are consistent with the given upper and lower bounds. 

§For Comet Halley the errors on the C/Si and N/Si ratios are dominated by our assumed factor of 2 error on the Si measurement. For 
the C/N ratio we assume a 50% error. 

^Nitrogen abundances relative to Si and C/N ratio provided by J. Raymond and given here. 
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